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Abstract 

The light fragment production reactions for 10 MeV to 3 GeV protons in- 
cident on 16 0, 27 A1, nat Fe, 93 Nb, and nat Ag are analyzed by a combination 
of an intranuclear cascade model and a generalized evaporation model which 
includes light nuclei up to Mg as ejectiles. It is concluded that evaporation 
is the dominant process by which particles lighter than or equal to Be are 
produced from targets heavier than O. 
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I. INTRODUCTION 



Fragment and residual nuclei production has attracted many people's interests, not only 
nuclear physicists but also astrophysicists and nuclear engineers. The production of in- 
termediate mass fragments from high energy proton-nucleus reactions or nucleus-nucleus 
reactions has been a hot topic in nuclear physics for a decade |l|J^] . Astrophysics and cosmic 
ray physics have been interested in residual nuclei production in order to calculate the pro- 
duction of cosmogenic nuclides in extraterrestrial matter by solar and galactic cosmic rays. 
Recently, nuclear engineering has needed the particle production cross sections for the devel- 
opment of accelerator-based systems for transmutation of radioactive nuclear waste. From 
the radiation safety aspect, it has also become more important to estimate the amount of 
radioactivity produced from various targets, as new applications of high energy proton ac- 
celerators, such as spallation neutron sources and the production of beams of unstable nuclei 
are being developed. 

In 1997, the Organization for Economic Cooperation and Development (OECD) Nuclear 
Energy Agency (NEA) conducted benchmark calculations on activation yields to determine 
the predictive power of current nuclear reaction models and codes [|J . The results calculated 
using many different codes which are based on a combination of different models, such as 
the intranuclear cascade model (INC), the exciton model, the evaporation-fission model, 
the quantum molecular dynamics model and the statistical multifragmentation model, were 
compared with experimental data. It became clear that most of the computer codes did not 
reproduce light fragment production reactions, such as Fe(p,X) 7 Be, especially at low proton- 
incident energy. They considered that an adequate description of the Fermi break-up model 
||] and a fragmentation model was urgently needed. 

The evaporation model has been very successful in describing residual nuclei production 
from hot nuclei. In many codes, not only those codes that were used in the benchmark cal- 
culation by OECD/NEA, but also the codes that are widely used for shielding calculations, 
such as the LAHET code ||, the evaporation model is used to describe the de-excitation 
of thermalized nuclei. Despite its success, the model has not been used to describe light 
fragment emission, except for a few studies concerning break-up of highly excited nuclei 

ML 

In this study, we propose a generalized evaporation model (GEVAP) for Monte Carlo 
simulation, based on the Weisskopf-Ewing model |7]||. Nucleons and helium nuclei are the 
dominant particles emitted from an excited nucleus. Therefore, only these particles are 
treated as ejectiles in the Dostrovsky's evaporation models implemented in the LAHET 
code H . On the other hand, some studies consider light nuclei heavier than a particles 
as ejectiles since there is no reason that those particles can not be emitted from excited 
nuclei via evaporation process. In our generalized evaporation model, 66 nuclides up to Mg 
are included as ejectiles, not only in their ground states but also in their excited states. 
Besides, we use the accurate level density function for the total decay width calculation 
instead of an approximate form of level density function which is used in the Dostrovsky's 
evaporation models 0. 

Light fragments produced from proton-induced reactions are analyzed by the combina- 
tion of the INC model implemented in the LAHET code H and the generalized evaporation 
model (GEVAP). In order to estimate light particle production from a nucleon-nucleus re- 
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action by the generalized evaporation model, we have to assume the ensemble of hot nuclei 
which are produced after the initial non-equilibrium stage. Since the excitation energy, 
the mass, and the charge of the hot residual nuclei produced from high energy reaction 
are widely distributed, we can not use the simple assumption that a single excited nucleus 
represents the ensemble of hot thermalized nuclei. However, the INC model can provide 
an ensemble of residual nuclei with broad distribution in excitation energy, nuclear mass 
and charge. The LAHET code employs the Bertini intranuclear cascade model [ID| for a 
non-equilibrium stage of nuclear reaction, and the Fermi break-up model [|J and the evap- 
oration model proposed by Dostrovsky et al. for a thermalized stage. Mass A±, charge 
Zj, excitation energy E, recoil energy, and the direction of recoil motion are extracted from 
the INC calculation done by the LAHET code. Then the de-excitation process of the hot 
nucleus with these quantities are calculated by GEVAP, instead of by the Fermi break-up 
model and the evaporation model employed in the LAHET code. In the following, we call 
this calculation procedure 'INC/GEVAP'. 

We focus mainly on 7 Be produced by proton-induced reactions in the energy range from 
10 MeV to 3 GeV, because 7 Be is the most intensively measured light fragment produced 
from various targets, and many experimental data are available for comparison. We compare 
the INC/GEVAP results with experimental data as well as the results calculated by using 
LAHET, to make the effect of using different de-excitation models clear. 



II. THE GENERALIZED EVAPORATION MODEL 

Let us consider that a parent nucleus i with an excited energy i£[MeV], a mass number 
Ai, and a charge number Zj emits a particle j in its ground state with Aj and Zj, and 
becomes a daughter nucleus d with Ad and Zd- According to the Weisskopf's formulation 
0, the decay probability Pj with total kinetic energy in the center-of-mass system between 
e and e + de is expressed as 

r> / \ j p d {E -Q-e) 

p A e ) de = 9j<7inv{£) tt^ ede, (1) 

Pi(E) 

where Oi nv is the cross section for the inverse reaction, pi and pd are level densities [MeV -1 ] 
of the parent and the daughter nucleus, respectively. With the spin Sj and the mass rrij of 
the emitted particle j, gj is expressed as gj = (2Sj + l)mj/TT 2 h 2 . In this study we use the 
Audi-Wapstra mass table to calculate the Q- values Q for emission of particle j. 



The cross section for the inverse reaction G{ nv is expressed as || 

a- (e) = { a9 ° n ^ + 6 / 6 ) for neutrons = a a (i _|_ d\ (2) 

m | UgCj (1 —V/e) for charged particles ~ 9 \ e J ' 

where a g = nR 2 [fm 2 ] is the geometric cross section, and V = ZjZ d e 2 /R c is the Coulomb 
barrier. 

In this study, we use the parameter set determined by Dostrovsky et al. || and Matsuse et 
al jr2"| . Dostrovsky et al. determined c n , Cj, b, Rb, and R c for n, p, d, t, 3 He, and a emission 
by fitting the expression to the theoretical calculation done by Shapiro Jl3[ and Blatt and 
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Weisskopf [[L4j], so that the effect of overlapping wave functions was taken into account. These 



parameters are used in the Dostrovsky's evaporation model || implemented in LAHET || . 
Meanwhile, Matsuse et al. determined the critical distance (Rf, and R c , with Cj = 1) by 
fitting Eq. (||]) to experimental fusion cross sections for heavy ion reactions. We use the 
Dostrovsky's parameters for n, p, d, t, 3 He, and a emission and the Matsuse's parameters 
for other particles. In the following we call these parameters "the precise parameter set". 
Besides the calculation with the precise parameter set, we use the simple parameter set, 
given by c n — Cj — 1, b — and Rf, = R c = r (A^ 3 + A]/ 3 ) [fm] for the inverse cross section. 
In the calculation with the simple parameter set, values of r = 1.2, 1.5, and 2.0 are tried 
to test the stability of our model. 

The total decay width Tj can be calculated by integrating Eq. flU) with respect to the 
total kinetic energy e from the Coulomb barrier V up to the maximum possible value (E—Q). 
By using Eq. (H) for a inv , the total decay width for the particle emission is expressed as 

r ' = ^f° e (l + f )"(*-«-<)*. (3) 

According to the Fermi-gas model, the total level density p(E) of a nucleus summed over 
all the possible states with the angular momenta is given by the expression [TH 



2y/a(E-S) 

where a = Ad/8 [MeV -1 ] is the level density parameter, and 5[MeV] is the pairing energy of 
the daughter nucleus evaluated by Cook et al. ||16|1 . For those values not evaluated by Cook 
et al, 5 obtained by Gilbert and Cameron |0J are used. E x is determined by Gilbert and 



Cameron JT5| as E x = U x + 5 where U x = 2.5 + 150/^4^. In the calculation with the precise 
parameter set, we use the Gilbert-Cameron-Cook-Ignatyuk (GCCI) level density parameter 
0, in which the pairing corrections and the energy dependence of the level density parameter 
are taken into account, instead of the simple expression a = A d /8. The GCCI level density 
parameter is employed in the LAHET code |J. 

When E is below E x , instead of Eq. (|4|) the following formula gives a good fit to the 



experimental level densities |15| : 



p{E) = L^~ E ^ T for E < E m (5) 



T 



where T is the nuclear temperature given by 1/T = ^Ja/U x — 1.5/U X . To connect Eq. 



and Eq. @ smoothly, E is defined as E = E x — T(log T — 0.25 log a — 1 .25 log U x + 2y/aU x ) . 
We use the expressions Eq. (§) and Eq. (13) to calculate the total decay width. The 



simple form p oc exp(2^ a(E — 5)), which is used in the Dostrovsky's evaporation models 
0, is a good approximation when the residual excitation energy is high, however, it is not 
applicable for residual nuclei with small mass and low excitation energy. 

When E — Q — V is below E x , Eq. (|3|) can be solved analytically, by substituting Eq. ([|) 
into Eq. (H). 
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r, = Zjr&ihit, t) + ((3 + V)I (t)} for E - Q - V < E x , (6) 

where Io(t) and I\(t,t x ) are expressed as: 

Io(t)=e-^ T (e*-l), 
h(t,t x ) = e- E °/ T T{{t - t x + l)e<* - t - 1)}, 

where t = (E — Q — V)/T and t x = E x /T. When E — Q — V is greater than E x , the integral 
of Eq. (|3p can not be solved analytically because of the denominator in Eq. ([|). However, 
it is expressed approximately as 

r i = TTTT^ **) + J 3(«, s*)e* + (0 + V) {I (t x ) + I 2 (s, s x )e s }} for E - Q - V > E x . 

(7) 

where I 2 (s, s x ) and J 3 (s, s x ) are given by: 



/ 2 (s, s x ) = 2V2 {s- 3 / 2 + 1.5s" 5 / 2 + 3.75s- 7 / 2 - (s; 3 / 2 + 1.5s; 5 / 2 + 3.75^ 7 / 2 )e s ^ s } , 
I 3 {s, s x ) = (V2a)- 1 [2s- 1 / 2 + As' 3 / 2 + 13.5 S - 5 / 2 + 60.0 S - 7 / 2 + 325.125 S - 9 / 2 - {(s 2 - s 2 x ) S - 
+ (1.5s 2 + 0.5s 2 )s; 5/2 + (3.75s 2 + 0.25s 2 )s; 7/2 + (12.875s 2 + 0.625s 2 )s; 9/2 
+ (59.0625s 2 + 0.9375s 2 )s^ 11/2 + (324.8s 2 + 3.28s 2 )s; 13/2 } e s *~ 
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with s = 2sJa(E - Q - V - 5) and s x = 2^Ja(E x - 5). 

In the present Monte Carlo simulation, ejectile j is selected according to the probability 
distribution calculated as pj = Tj/^jFj, where Tj is given by Eqs. @ or (|^). The total 
kinetic energy e of the emitted particle j and the daughter nucleus is chosen according to 
the probability distribution given by Eq. (JJ). The angular distribution of the motion is 
randomly selected to be isotropic in the center-of-mass system. The excitation energy of the 
daughter nucleus Ed is calculated as Ed = E — Q — e. 

In this study, we consider 66 nuclides as ejectiles, not only in their ground states but also 
in their excited states. It is important to include excited states in the particles emitted via 
the evaporation process, because it greatly enhances the yield of heavy particle emission . 
The selected ejectiles satisfy the following criteria: (1) isotopes with Zj < 12; (2) naturally 
existing isotopes or isotopes near the stability line; (3) isotopes with half-life longer than 1 
ms. The selected ejectiles are listed in Table |. 

If the mean lifetime of a resonance is longer than the decay width of the resonance 
emission, such a resonance can survive during the evaporation process. The excited state is 
included if its half lifetime Ti/ 2 [sec] satisfies the following condition: 

Tl/2 -> h 

to2 > ' (8) 

where T* is the decay width of the resonance emission. T* can be calculated in the same 
manner as for a ground state particle emission. The Q-value for the resonance emission is 
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expressed as Q* = Q + E*, where E* is the excitation energy of the resonance. The spin 
state of the resonance S* is used in the calculation of gj, instead of the spin of the ground 
state Sj. We use the ground state mass rrij for excited states because the difference between 
the masses is negligible. 

Instead of treating a resonance as an independent particle, we simply enhance the decay 
width of the ground state particle emission. We redefine the decay width Tj as 

ri = r? + £r» (9) 

n 

where T° is the decay width of the ground state particle j emission, and is that of the 
nth excited state of the particle j emission which satisfies Eq. (||). 

The total kinetic energy distribution of the excited particle emission is assumed to be the 
same as that of the ground state particle emission. Sj, E*, and Ti/ 2 used in this study are 
extracted from the Evaluated Nuclear Structure Data File (ENSDF) database maintained 



by the National Nuclear Data Center [17 



III. COMPARISON WITH EXPERIMENTAL DATA 

The excitation functions of 7 Be produced by proton reactions on 16 0, 27 A1, nat Fe, and 
93 Nb are shown in Fig. |l| The results calculated by INC/GEVAP with the precise parameter 
set, which consists of the parameters for inverse reactions determined by Dostrovsky et al. 
and Matsuse et al |jl2[ and the GCCI level density parameter, are shown by the solid lines. 
The estimates by INC/GEVAP with r$ = 1.5 and those by LAHET are also shown in the 
figures as well as the experimental data collected in Ref. 0. The results by INC/GEVAP 
with r = 1.2 and 2.0 are represented only for the Nb target, because there is less difference 
in the results for other targets. For the O target, the differences in the estimates between 
by INC/GEVAP with r = 1.5 and by that with r = 1.2 or 1.5 are 20 %, except at the 
threshold energy. For Al target, the differences are within a factor of two in the whole 
energy region. INC/GEVAP with r$ = 2.0 produces almost the same cross sections for the 
Fe(p,X) 7 Be reaction as those with precise parameter set, and the differences are within 20 
%. In the whole energy region, for all the targets, INC/GEVAP produces more 7 Be as r 
increases. The estimates by INC/GEVAP for Al with r = 1.5 give the best agreement with 
the experimental data, as seen the dashed line lying underneath the measurement points 
in the whole energy region in Fig. [I] (b). Whereas for Fe and Nb INC/GEVAP with the 
precise parameter set reproduce the excitation functions better than that with the simple 
parameter set, and the estimates agree with most of the experimental data within 50 %. 

INC/GEVAP reproduces the excitation functions for all the targets, whereas LAHET 
fails to reproduce the shape of the excitation functions except for the O target. The shapes 
of the excitation functions estimated by INC/GEVAP do not change with the choice of the 
parameter sets. Since LAHET severely underestimates the 7 Be productions from Al below 
300 MeV, Fe and Nb below 3 GeV, it is obvious that the Fermi break-up is not the dominant 
process for the 7 Be productions in these reactions. 

The isotopic distributions of H, He, Li, and Be nuclei produced from 2100 MeV proton 
incident on 16 and from 480 MeV proton incident on nat Ag are shown in Fig [2|. The results 
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estimated by INC/GEVAP with the precise parameter set (the open squares) are shown 
as well as the experimental data for the O target measured by Olsen et al. | pL 8|| , and the 
data for the Ag target by Green et al. [[19] (the closed circles). INC/GEVAP reproduces 
the isotopic distributions for both these reactions, and the estimates agree with most of the 
measurements with 50 % accuracy 

IV. SUMMARY AND CONCLUSION 

We have formulated a generalized evaporation model (GEVAP) based on the Weisskopf- 
Ewing model |7]||. The features of the model are: (1) the accurate level density function 
is used for deriving the decay width of particle emission; (2) sixty-six nuclides up to Mg, 
not only in their ground state but also in their excited states, are taken into account in this 
study 

The combination of the intranuclear cascade model (INC) and GEVAP successfully re- 
produces the excitation functions of 7 Be produced by protons incident on 16 0, 27 A1, nat Fe, 
and 93 Nb. The choice of the parameter set in GEVAP does not affect the resulting shapes 
of the excitation functions. INC/GEVAP also predicts the isotopic distributions of H, He, 
Li, and Be produced from O and Ag with 50 % accuracy. 

From the results, it is concluded that the evaporation process is the main process via 
which particles lighter than or equal to Be are produced by protons incident on targets 
heavier than O. INC/GEVAP can predict the production cross sections of particles lighter 
or equal to Be between 50 % to a factor of two in accuracy depending on the choice of the 
parameter set used in GEVAP. In this study, the precise parameter set gives the best results 
for overall reactions, and the accuracy is 50 % on the average. 
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Proton incident energy (MeV) Proton incident energy (MeV) 



FIG. 1. The excitation function of 7 Be produced from 16 0, 27 Al, nat Fe and 93 Nb: The estimates 
by INC/GEVAP with ro = 1.5 and the precise parameter set are shown by the dashed lines and 
the solid lines. The results calculated by using LAHET are shown by the dotted lines. The open 
squares are the experimental data collected in Ref. ||], and the crosses and the open diamonds are 



experimental data in Ref. pQ] and Ref. [21], respectively 



(a) 0(p,X), E = 2100 MeV 



(b) Ag(p,X), E = 480 MeV 




2 3 3 4 5 
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FIG. 2. The isotopic distributions of the nuclei produced from: (a) 2100 MeV protons incident 
on 16 0; (b) 480 MeV protons incident on nat Ag. The open squares with the solid lines denotes the 
results by calculated INC/GEVAP with the precise parameter set, and the black circles with the 
dashed lines denotes experimental data [ 18, 1^] . The lines are drawn to guide to eyes. 
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TABLES 

TABLE I. The ejectiles considered in this study 



Zj Ejectiles 






n 














1 


p 


d 


t 










2 


3 He 


4 He 


6 He 


8 He 








3 


6 Li 


7 Li 


8 Li 


9 Li 








4 


7 Be 


9 Be 


10 Be 


n Be 


12 Be 






5 


8 B 


io B 


n B 


12 B 


13 B 






6 


io c 


n c 


12 C 


13 C 


14 C 


15 C 


16 C 


7 


12 N 


13 N 


14 N 


15 N 


16 N 


17 N 




8 


14Q 


15Q 


16Q 


17 


18Q 


19Q 


20Q 


9 


17 F 


18p 


19p 




21p 






10 


18 Ne 


19 Ne 


20 Ne 


21 Ne 


22 Ne 


13 Ne 


24 Ne 


11 


21 Na 


22 Na 


23 Na 


24 Na 


25 Na 






12 


22 Mg 


23 Mg 


24 Mg 


25 Mg 


26 Mg 


27 Mg 


28 Mg 



10 



